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We consider dynamics of classical dimer models undergoing a phase transition to an ordered, 
frozen state. Relaxation processes are dominated by traps which are entropic in origin and can be 
traced to the locally jammed nature of the dimer states. Depending on the nature of the phase 
transition, critical dynamics are characterized either by an exponential, or sub-exponential in time 
relaxation of the order parameter. In the latter case relaxation time scales diverge exponentially as 
the critical point is approached, following the Vogel-Fulcher law. 
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A variety of systems such as supercooled liquids, col- 
loids, granular matter and foams, exhibit a transition 
from a flowing fluid phase to a frozen solid phase. Jam- 
ming due to spatial constraints imposed on the elemen- 
tary constituents of these materials has been proposed as 
a possible common cause of this dynamical arrest [Q-^. 
Model systems, such as hard spheres, have an impor- 
tant role to play in the investigation of such a scenario 
since they allow for a precise definition of jamming [Q. 
They are also useful in elucidating the precise relation- 
ship between thermodynamics and dynamics in materials 
exhibiting a jammed phase |5j. 



Dimer models are examples of jammed systems which 
have the added advantage of being exactly solvable 
States of the dimer model are specified by placing dimers 
on the bonds of the lattice so that every site is covered by 
exactly one dimer; see Fig |^. These dimer coverings are 
"locally jammed" Q as each dimer cannot move to an 
empty neighboring bond without violating the packing 
constraint. Moves that involve loops of dimers and adja- 
cent empty bonds, on the other hand, arc allowed. An ex- 
ample of such a move for the hexagonal lattice involving 
an elementary plaquette is shown in Fig. (Stochastic 
dynamics of the dimer model on the square lattice based 
on these elementary moves were considered by Henley 
0.) Most states of the dimer model allow for elemen- 
tary moves; an example of one which does not is shown 
in Fig. The smallest change in this case involves a 
system spanning loop, and we call this state "maximally 
jammed" . If we define an energy functional on the space 
of dimer coverings which favors the maximally jammed 
state, a transition into this state will occur as we lower 
the temperature. The central question we address in this 
letter is: What happens to relaxation time scales of the 
dimer model as the transition to the maximally jammed 
state is approached? We will show that the relaxation 
is dominated by entropy barriers and is sensitive to the 
thermodynamics of the model near the phase transition 
point. 





FIG. 1. (a) Dimer covering of the honeycomb lattice with 
an elementary loop update indicated by the arrows. The 
numbers are the heights of the equivalent interface, (b) An 
ordered, maximally jammed dimer covering; the equivalent 
interface is tilted with maximum slope. 

We consider two energy functionals. One leads to 
a Pokrovsky-Talapov type transition to the maximally 
jammed state ||^. The other exhibits a continuous tran- 
sition to the same state along a metastable line, and be- 
longs to a different universality class. For both tran- 
sitions we find strong departures from the canonical 
critical-slowing-down scenarios |^ due to the presence of 
entropy barriers. Barriers can be traced directly to the 
non-local nature of the dynamical moves allowed by the 
jammed states. For the transition that occurs along a 
metastable line, time scales associated with critical fluc- 
tuations diverge exponentially following a Vogel-Fulcher- 
like form. This is reminiscent of what is observed in frag- 
ile glass formers 1^. To our knowledge this is the first 
example of a system with no quenched in disorder which 
exhibits this type of critical dynamics. 

a. Dimer models We consider the dimer model on 
the 2-d hexagonal lattice of linear size L, having 2L^ 
sites and bonds, with periodic boundary conditions 
Ip^ . A useful representation of the dimer model is given 
by the height map which associates a discrete interface 
h{x, y) with every dimer covering 0. The heights of the 
interface are defined on the vertices of the dual triangu- 
lar lattice. The height difference A between two nearest 
neighboring sites is -2 or -1-1 depending on whether the 
bond of the honeycomb lattice that separates them is oc- 
cupied by a dimer or not; see Fig. |^a. Directions in which 
the height change is -l-A are specified by orienting all the 
up pointing triangles of the dual lattice clockwise. 

The dimer model has an extensive entropy. The ensem- 
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ble of equal weighted dimer coverings maps to a rough 
surface with a gradient-square free energy 0. Fluctua- 
tions of the surface are entropic in origin. A phase tran- 
sition can be induced in the dimer model by including an 
energy functional which is minimized by a dimer cover- 
ing corresponding to a smooth, maximally tilted surface 
which corresponds to the maximally jammed state shown 
m Fig. |l|b. Here we discuss two such energy functionals. 

For periodic boundary conditions the tilt vector, 
{Axh, Ayh), where A^^yh is the average height difference 
in the x or y direction, has only one independent compo- 
nent p . In terms of p, the energy functionals can be 
written as: 




fi is a dimensionless coupling that drives the transition. 

The entropy of the dimer model as a function of p was 
calculated exactly pO| , [T3t : 

^-(l-p) + -y dxln[cos:E]^ (2) 

This function has a maximum at p = which is the 
equilibrium value a,t n = 0. For finite fi the / model was 
previously considered in Ref. , while the K model was 
solved by Kasteleyn |Q . 

In the / model, with the free energy F = E — S, and 
the energy and entropy given by Eqs. |^ and |^, there is 
an interesting phase transition along the metastable line, 
when the order parameter is confined to the free energy 
well around the zero-tilt state. Namely, at p^, = it/{8^/3) 
(the end-point of the metastable line) the order parame- 
ter p has a discontinuous jump from to 1, characteristic 
of a first-order transition. At the same time, as /i* is ap- 
proached from below, fiuctuations of p around diverge, 
as would be expected at a critical point. This transition 
was discussed in detail in Ref. 

By contrast, in the K model, there is a continuous 
variation of the equilibrium value of the order parame- 
ter, Pmim with fi. On varying p from to the transition 
point p^ = ln2, Pmin varies from to 1. Apart from 
the shift of the free energy minimum, the fiuctuations 
also grow with increasing p, and at the transition, the 
curvature of the free energy at its minimum goes to zero. 

b. Coarse-grained dynamics As mentioned in the in- 
troduction, the hard constraint of non-overlapping of 
dimers, gives rise to nonlocal dynamics. We consider 
Monte-Carlo dynamics based on loop updates with loops 
of arbitrary size; a concrete implementation is given in 
Ref. [ p^ . Since we take periodic boundary conditions, 
loops with different winding numbers can be formed. 
We restrict loop updates to loops with winding numbers 
(0,0), (1,0) and (0,1), only. The microscopic transition 
rates for loop updates are given by Metropolis rules that 
follow from Eq. n. 



Given the microscopic loop dynamics, which satisfy 
conditions of ergodicity and detailed balance, we ask 
what are the coarse-grained dynamics of the order pa- 
rameter, p. Since the energy functions in Eq. |l| depend 
on the global tilt p only, it follows that all updates of 
topologically trivial loops (i.e. those with (0, 0) winding 
number) have AE = 0. Only when system spanning 
loops with nonzero winding numbers are updated does 
the energy of the state change |l^ . This feature natu- 
rally leads to fast and slow processes in the Monte-Carlo 
dynamics. On a faster time scale, non- winding loops are 
updated with no effect on the overall tilt of the surface, 
while on a much slower time scale, winding loops are up- 
dated causing a change in the tilt of the surface. 

The coarse-grained dynamics of global tilt changes are 
described by a master equation, which features an un- 
usual form for the transition rates between different tilt 
states. Namely, the rates of transitions from higher into 
lower tilt states are determined by the energy barrier 
alone: 



here Fq is a constant. This follows from the observation 
that in order to lower the tilt and increase the energy an 
existing system spanning loop needs to be updated. 

From detailed balance, M^p-i/L,p/W^p,p-i/L — 
exp[— (F(p — 1/L) — F{p))], we conclude that the rates 
of transitions to higher tilt states must be determined by 
the entropy: 



Wp,p_i/i=r„e-(«(''-i/^)-S(rt) . (4) 



Eq. y can also be argued from the following observa- 
tion: to increase the tilt and lower the energy, a new 
system-spanning loop has to be accommodated by (typ- 
ically) many rearrangements of the topologically trivial 
loops. Once a new system spanning loop is introduced, 
and a higher tilt state is reached, the entropy decreases. 
The form of the transition rates that we are arguing for 
here, was directly observed in numerical simulations of 
the three coloring model |l7j , which is a close relative of 
the dimer model. The two are equivalent if, in the dimer 
models, a weight of 2 is attached to each loop formed by 
bonds that are not covered by dimers. 
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mean that once the system decays out of state p (in time 
Tp) it quickly looses memory of where it came from. 

To quantify the tilt dynamics wc compute the tilt-tilt 
autocorrelation function C(i), defined as: 



C(t) 



{piof) - {piO)f 



(6) 




FIG. 2. (a) The time scales for relaxing out of different tilt 
states p for the I and K models (scaled by L), for values of p 
below the transition, (b) The tilt-tilt autocorrelation function 
for the two models. The full line is obtained from Eq. ^ while 
the dashed line is a result of the saddle point evaluation of 
Eq. ^. L was chosen to be 4096 (/) and 24 (K) to make the 
time scales comparable. Time is measured in units of ^. 

The first consequence of the above form of the tran- 
sition rates is that the time scale of relaxation out of 



p+l/L,p 



a state with tilt p, Tp — l/(W^p-i/i,p 
non- decreasing function of p. The exact expressions for 
Tp (measured in units of T"^), 



-i§-ppL I -L[|ln2+f ln[cos(f -^)]] 
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g-fMi + g-L[fln2 + §ln[cos(f -fi)]] 



(I) 
(K) 



(5) 



follow from Eqs. || and ^, and they are plotted in Fig. ||. 
For the / model the time scale increases monotonically 
with p as in the hierarchical models of Palmer et 
al. while in the K model an initial rise is followed 
by saturation beyond Pmm- Both are in sharp contrast 
with Langevin dynamics around the equilibrium state, 
for which the time to relax out of a given macro-state 
decreases the further the order parameter is away from 
its equilibrium value. (For example, in the Ising model, 
with Glauber dynamics and in the disordered phase, the 
relaxation time out of a given magnetization state de- 
creases with increasing magnetization.) 

The exponential separation between different Tp's for 
all p in the / model, and for p < pmin in the K model, 
implies that the p dynamics are trap-like pU . By this we 



with the average taken over different histories of p. As- 
suming trap-like dynamics. 



Cit) 



'^p-F{p)p-t/T^ 



(7) 



i.e., C{t) is an equilibrium weighted average of relax- 
ations out of different p states. Such an approximate 
model works very well when checked against numerical 
simulations of the three-coloring model |17|. 



The asymptotic decay of the autocorrelation functions 
obtained from Eq. ^ can be extracted from a saddle point 
analysis of the sum and using a quadratic approximation 
for the entropy (Eq. ||). The autocorrelation functions 
obtained from the sum are compared to the saddle point 
solutions in Fig. ||b. The latter can be used to analyze 
the asymptotic form of C{t), and the dependence of the 
relaxation time scales on p. In the limit of t —^ oo, and 
up to logarithmic corrections in t, we find: 



C{t) 



exp{- 



{I) 



^[ln(^)]n {K) 



(8) 



showing that C{t) has an exponential decay in the K 
model, and a slower than exponential decay in the / 
model. From Eq. || we also conclude that the relaxation 
times scale, r, for the decay of C{t) to an arbitrary con- 
stant Co, diverges exponentially as p ^ p<t for the / 
model. This is the Vogel-Fulcher type behavior observed 
in many fragile glass formers. First order corrections to 
Eq. ^ lead to an even more rapid increase of time scales, 
with r/lnr diverging as Vogel-Fulcher. In the K model 
the corresponding time scales are exponential in p (Ar- 
rhenius form) with no divergence at the transition point. 
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FIG. 3. Barrier height, B{p) (dimensionless) shown as a set 
of solid hnes, and the quadratic approximation to the dimen- 
sionless free energies of the / (dotted line) and K (dashed 
line) model; n is chosen close to /i, and L — 24. Note the 
logarithmic scale for the barrier height. 



these barriers (see Fig. ||). 

In the / model the free-energy minimum is centered 
at p = for all values of /i. The well gets shallower 
as fi^ is approached implying that the system explores 
larger and larger barriers leading to a divergent relax- 
ation time scale. In the K model, the minimum itself 
shifts to larger values of p and regions of increasing bar- 
rier heights {pmin — 1 in Fig. ^). The dynamics is then 
completely dominated by a single time scale determined 
by the barrier at Pmin- 

In conclusion, we have shown that classical dimer mod- 
els with loop dynamics are an interesting model system 
for studying glass-like dynamical arrest close to a critical 
point. Our analysis shows that the jammed nature of the 
dimer states plays a crucial role in determining the nature 
of relaxation processes. To quantify how this is reflected 
in the relaxation time scales we introduced an effective 
dynamical model based on the master equation for the 
order parameter. We expect that similar ideas can be 
carried over to the investigations of hard sphere systems 
which are of direct relevance to colloidal experiments . 
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The coarse grained dynamics defined by the transi- 
tion matrix elements, Eqs. || and |[ were argued to fol- 
low from the nonlocal loop dynamics of the dimer mod- 
els. From this form of the W^-matrix all the conclusions 
about critical dynamics of the / and K model are de- 
rived. We have confirmed this picture in considerable 
detail in simulations of the three coloring model |^,^ , 
which, as discussed earlier, is the loop weighted dimer 
model. The loop weights are not expected to affect the 
qualitative features of the energy and entropy function- 
als. Indeed, the measured Tp for the / and K cases of 
the three-coloring model compare very well |2^j2^ to the 
analytical form plotted in Fig. |2| Similarly, the tilt-tilt 
correlations show exponential decay in the K model ||2^ ] 
and sub-exponential decay in the / model [21 2^ as ex- 
pected from the analysis presented in this paper. The 
numerical evidence for Vogel-Fulcher type divergence of 
the relaxation time scale in the / model was reported 
previously [2l| ]. 

The origin of the difference in the dynamical behavior 
of the / and K models can be traced back to the inter- 
play between the free energy and the dynamical barriers. 
The transition rates presented in Eqs. ^ and 0, can be 
interpreted in terms of a barrier B{p) ~ e('S{p-i7i)--S(p)) 
dividing the usual Metropolis rates defined in terms of 
the free energy and leading to Langevin dynamics [ p^ . 
These barriers increase exponentially with p as illustrated 
in Fig. ^. Dynamics of the order parameter can be viewed 
as relaxation in the free energy well in the presence of 
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